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Abstract. Relativistic thermodynamics is treated from the point of view of kinetic theory. It is shown 
that the generalized Jiittner distribution suggested in [T] is compatible with kinetic equilibrium. The 
requirement of compatibility of kinetic and thermodynamic equilibrium reveals several generalizations 
of the Gibbs relation where the velocity field is an independent thermodynamic variable. 



1 Introduction 

Relativistic thermodynamics is the less established a- 
mong relativistic theories. The famous "imbroglio" of 
the temperature of moving bodies |2] indicate, that 
basic thermodynamic concepts may have contradic- 
tory explanation in relativity theory. The first general- 
ization of thermodynamics in a relativistic framework 
was the proposal of Planck and Einstein [2111], and 
were questioned by several critical approaches more 
than half century later (SJIHIIZIIH] ■ These suggestions 
seemingly exclude each other, but their arguments 
are convincing enough to provoke a longstanding dis- 
cussion that continues even today [51 [TU1[TT1[T^I131II41 
[TSllTCllTTlfTSlfT^ . On the other hand the experimental 
measurement of temperature of fast moving bodies 
is an everyday practice in astrophysics and cosmol- 
ogy PPII^ and in heavy ion physics [^[^ and one 
can deal with thermal aspects seemingly without ap- 
parent problems. From mentioned investigations it is 
clear that the problem requires more than a covariant 
form of the physical quantities and the Gibbs rela- 
tion, or a definition of thermodynamic bodies (mov- 
ing boundaries, etc.). The interpretation of thermody- 
namic equilibrium among interacting bodies seems to 
be the most important aspect. 

The history of temperature of moving bodies shows, 
that the problem cannot be analysed by argumenta- 
tion based exclusively on thermodynamics and spe- 
cial relativity. Moreover, statistical mechanics alone 
does not give a definite answer, statistical explana- 
tions were given to all of the previous different views 
providing different interpretations of the macroscopic 
quantities P^[5]. The reason is that the problematic 
point is somehow outside the usual realm of equi- 
librium theories. The key aspect is the interaction 
of moving bodies, the thermodynamics of the mo- 
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tion. Therefore the problem is best investigated in the 
framework of classical field theories - like continuum 
(hydrodynamics) or kinetic theories - where this as- 
pect is the part of the theory and equilibrium and 
thermodynamic equilibrium is motion related. How- 
ever, one should be aware, that both in hydrodynam- 
ics or in kinetic theory the definition of equilibrium is 
a starting point where several basic physical quantities 
are interpreted, and thermodynamics and thermody- 
namic equilibrium plays a key role. 

In our previous research we investigated stabil- 
ity and causality in relativistic hydrodynamics. A de- 
tailed analysis of the Second Law revealed [24: that as- 
suming acceleration independent entropy production 
the local rest frame entropy density may depend on 
the absolute value of the energy momentum vector 
s{^/WE^,n), where = UbT°-^ , with m° is the four 
velocity of the fluid and T"** is the energy-momentum 
density tensor. Then the consequent Gibbs relation 
can be written as 

de + — dg" Tds + iidn. (1) 

e 

Here e — UaT'^^'ui, is the local rest frame energy 
density, and Qa — AabE^ is the momentum density. In 
the usual treatment the second term of the left hand 
side is missing, the entropy depends only on the en- 
ergy density and particle number density s(e, n). Our 
generalization eliminated the generic instability of the 
original Eckart theory [25] and we have proved that 
the stability is independent of any flow frames like 
Eckart or Landau-Lifshitz. The conditions are the non- 
negativity of the transport coefficients and thermody- 
namic stability, the concavity of the entropy density. 
This is in strong contrast to the usual theory of Eckart, 
which is unstable [5^, or to the Israel- Stewart theory, 
where there are several complicated and counterintu- 
itive conditions beyond the above mentioned natural 
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requirements [271128) . Other different recent investiga- 
tions do not give similar stability properties [2911301 
[31]. 

Motivated by the previous results, and with the 
help of a generalized internal energy and modified ba- 
sic thermodynamic relations, we have analysed the 
problem of relativistic temperature transformations 
from the point of view of relativistic hydrodynamics. 
A homogenization of the energy-momentum balance 
shows that the entropy can depend directly on the 
energy- momentum vector of the body S{E'^,V) and 
results in the following form of the Gibbs-relation: 



gadE"" = TdS - pdV, 



(2) 



where S and V are the entropy and volume of the 
body and is the four- vector intensive quantity con- 
jugated to the total energy-momentum four vector. 
Here we did not consider changes in particle number 
and assumed a constant velocity. The thermodynamic 
body was interpreted as a homogeneous continuum, 
where the local rest frame is synchronized and we re- 
quired also a constant average velocity. The above re- 
lation implies that the entropy may depend on the en- 
ergy momentum vector and not only on its absolute 
value. We have deduced that g° is timelike and not 
necessarily parallel to the velocity of the body. The 
appearance of a perpendicular component resulted in 
nontrivial conditions for thermodynamic equilibrium 
of two thermodynamic bodies. In case of a single body 
one may sit into a preferred frame, determined by 
the energy-momentum vector, but for interacting bod- 
ies we should consider the energy-momentum of each, 
therefore we do not have this possibility. The condi- 
tion of equilibrium requires the equality of the tem- 
peratures and the composite velocities, determined by 
the vectorial intensive quantity. These velocities are 
not identical to the velocities of the respective bodies 
[1]. With this result we could give a physical back- 
ground of the historical temperature transformation 
formulas of Einstein-Planck, Blanusa-Ott and Lands- 
berg. Moreover, the experimentally observed Doppler- 
like transformation of the spectral temperatures was 
derived, too [3a[ia[l4]. 

The point of view of a related, but more general 
theory was essential. On the other hand the second 
form of the Gibbs relation © is more general than 
the first one ([T]): the entropy depends on the energy- 
momentum vector, not only on its absolute value. The- 
refore in this paper we look for an independent, but 
related theory, and investigate whether kinetic theory 
is compatible with the generalized thermodynamic re- 
lations HI) and/or In these researches we focus on 
the following aspects of the basic problems: What is 
the meaning of the rest frame in the form of the ther- 
modynamic relations? How can one distinguish dissi- 
pative and nondissipative parts of physical quantities? 
These questions appear as the choice of flow-frames 
(like Eckart or Landau-Lifshitz) in relativistic hydro- 
dynamics, arise as the existence of Brenner diffusion 



velocity in non-relativistic hydrodynamics |33) . are at 
the root of the temperature transformation paradox 
and in its most general form they are related to the 
very meaning of covariance and the existence of abso- 
lute objects in relativistic theories [M1I35] . 



2 Basic hydrodynamics and the separation 
of dissipative and non-dissipative 

In the following the Lorentz form is diag{l, —I, —I, —I) 
and the units are chosen that speed of the light and 
the Boltzmann constant are one c — 1, kB = 1- For 
one component fluids, the basic macroscopic fields are 
the particle number density four vector A^°, the en- 
tropy density four vector 5*° and the energy momen- 
tum density tensor T"''. The particle number and the 
energy-momentum are conserved in the following in- 
vestigations, therefore 



0, 
0. 



(3) 
(4) 



However, the entropy is not conserved, the en- 
tropy production is zero only in thermodynamic equi- 
librium: 

daS'^ix") > 0. (5) 

With the help of a velocity field u"'{x^) the particle 
number and entropy density four vectors and the en- 
ergy momentum density tensor can be expressed by 
local rest frame quantities as 





= nu" + 




= su" -1- 


rpab 


= E'^u'' 




= eu^v!' 







J", 



(6) 
(7) 

(8) 



Here the four component physical quantities are 
splitted to parallel and perpendicular components to 
the local velocity, timelike and spacelike components 
in the local rest frame, n :— N°-Ua is the particle 
number density, j° := A°-^Nb is the diffusion current, 
s := S°'Ua is the entropy density, J° := A°'^Nb is the 
entropy current, e :— UaT°'^Ub is the energy density, 
g" := A'l:T'*Ub is the momentum density (and the en- 
ergy current) and P'^^ := Z\JJZ\^T^'' is the pressure ten- 
sor. Here we have introduced the u-orthogonal projec- 
tion Z\°^ := 5""^ — The energy-momentum tensor 
is symmetric and therefore the pressure tensor is also 
symmetric and the momentum density and the en- 
ergy current are equal. Moreover, in addition to these 
usual local rest frame quantities we have introduced 
the energy-momentum vector E"" := UbT"^^ and the 
energy-momentum current density Q""^ :— A^T"^^ in 
order to clearly distinguish between momentum den- 
sity g" := A'^E^ and energy current := u^Q'^^. The 
forms ©-(El) conveniently express the particle number 
density four vector and the energy-momentum density 
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tensor with the help of local rest frame quantities rel- 
ative to a velocity field. 

The balances ©-([l]) can be written as 

daN'' = n + ndau'' + daf, (9) 

= eu" + eu^'dbu" + 9° + q^dbU^ + eu° + 
u'^dbq'' + q^dbu'' + dbP''\ (10) 

where the dot denotes the derivative along the world 
lines of the velocity field (substantial time derivative) 

It is considered as an evident fact that we can ad- 
ditively decompose the particle number density four 
vector and the energy-momentum tensor into a non- 
dissipative and a dissipative part. We usually assume, 
that these parts are easily distinguished in the local 
rest frame, because there the non-dissipative particle 
number density four vector is parallel to the veloc- 
ity and the non-dissipative energy momentum tensor 
Tq'' is diagonal: 



K 

-tab 



TIqU 



T^" = eouV - poA 



ab 



(11) 

(12) 



Here po is the static, scalar pressure, determined 
by the static equation of state of the fluid. Therefore 
the diffusion current density j'^ is the dissipative part 
of the particle number density and the momentum 
density /energy current and the difference of the 
total and the equilibrium pressure 77°'' — P'^'' +poA°-^ 
- the viscous pressure - are the dissipative parts of the 
energy momentum. 

Let us observe that with this splitting the distinc- 
tion of dissipative and non-dissipative parts of the 
physical quantities is related and depends on the lo- 
cal rest frame, the velocity field of the continuum and 
requires a particular thermostatics to determine the 
static pressure. However, the dissipation is principally 
defined by the entropy production, and therefore im- 
plicitly related to the background thermostatics, that 
is to the concept of local equilibrium. On the other 
hand, in case of dissipation the rest frame is not deter- 
mined by the local thermodynamic equilibrium, nei- 
ther by any special form of the physical quantities. 
One can get restrictions of the possible equilibrium 
entropy four vector considering simply that general 
information |36| . Moreover, one may consider to ex- 
tend the thermodynamic state space by the velocity 
field. In this case a pure continuum approach is based 
on the entropy inequality and requires some special 
mathematical methods (Coleman-Noll or Liu proce- 
dures [371I3H])- 

In this paper we choose an other way and con- 
sider kinetic theory as a general background. We will 
explore whether the previously mentioned generalized 
relativistic thermostatics introduced by pure phenome- 
nological arguments is compatible with kinetic theory. 



3 Equilibrium in kinetic theory 

For the sake of simplicity we consider here a single 
component Boltzmann gas without external fields, 
where the momentum vector of a particle is p° and its 
mass is m = y/p°-pa ■ The Boltzmann equation for the 
one particle distribution function f{x°',p°'), is 



P'daf = C{f). 



(13) 



Here da denotes the partial derivative of / regard- 
ing the spacetime event cc" and C(/) is the following 
collision integral: 

C{x,p) duj,du,'du[ {ff[W{p',p[\p,Pi)- 

ffiW{p,pi\p',p[)), (14) 

where duj = is the momentum space measure. 
W{p,pi\p' ,p[) denotes the transition rate from a par- 
ticle pair momentum {p'°',Pi) before the collision, to 
particle pair momentum {p'^,p\) after the collision. 

The particle number density vector, the energy- 
momentum density tensor and the entropy density 
vector are defined as 



(15) 
(16) 
(17) 



:= J dup'^f, 
T^" := J dLop^f, 

:= I dojp-f{\nf-l) 



The balances of these fields are calculated from the 
Boltzmann equation (see e.g. [39]), exploiting the prop- 
erty that due to the microscopic conservation laws a 
linear combination of a constant and the momentum 
four vector 'ip{x°',p°') = a{x"') + ba{x'^)p"' is a collision 
invariant, does not change in collisions: 



dw^C = 0. 



As a consequence of this property the conservation of 
the number of particles ^ and energy-momentum Q 
follow. 

We may calculate the entropy production, too. As- 
suming momentum independent force fields and uni- 
tary scattering matrices we get a local form of the 
H-theorem: 



das'" ^-j dulnfp^daf^ 

i J dududuidw[ - In ^ - 1^ > 0. (18) 

This entropy production vanishes if and only if the 
distribution function satisfies the simple functional re- 
lation: 

/(x^p"^)/(a;^pr) = /(a;^p'^)/(x^p?). (19) 
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We assume that the distribution function tends 
to a definite hmit as time progresses and the system 
develops into an equihbrium state where the entropy 
attains its maximum vahie. This condition, together 
with the requirement, that the equihbrium distribu- 
tion function must be a solution of the transport equa- 
tion uniquely determines its form. 

In case of binary collisions the momenta satisfies 
the conservation of energy-momentum 



P +Pl=P 



■Pi- 



rn 



Moreover, the local equilibrium distribution function 
/o satisfies the functional relation ([T9|) . Therefore the 
logarithm of the equilibrium distribution function is a 
collision invariant, and necessarily has the form 



(21) 



with arbitrary space- and time-dependent parameters 
a{x'') and /3°'{x^). Let us remark, that the equilibrium 
distribution function is not a solution of the transport 
equation without any further ado. Substituting into 
(|T5)) we get, that a and /?" must satisfy 



(22) 



Global equilibrium is defined by the requirement, that 
this equation must be an identity for arbitrary values 
of This implies simple initial (boundary) condi- 
tions and gives that a is a constant field and /3° obeys 
the Killing equation: 



a{x) = const, daP^ + dbP" = 0. 



(23) 



Local equilibrium considers more involved initial (bo- 
undary) conditions of the transport equation and the 
consequent macroscopic balances are used to interpret 
and evaluate the equilibrium fields through the equa- 
tions p^-pH). 

With the help of the equilibrium distribution func- 
tion we may calculate the particle number balance: 

= I dujp'^dafo - J diOp'^p'/Mb = 

^Ngdaa^T^'dafSb. (24) 

In order to get thermodynamic relations one can 
calculate the Boltzmann-Gibbs entropy with the equi- 
librium distribution (PTI) . too: 

SS^ J dc^pVo(ln/o-l) = (l-a)iVo" + W^ (25) 

(124)) and ^ results in 

daSS + adaK - f^bdaT"' = 0. (26) 

The above relations indicate a kind of Legendre 
transformation between the formulas (l24l) and (l26l). 



However, (j26p cannot be interpreted as a generalized 
Gibbs-relation, we cannot consider the divergences as 
differentials without any further ado. This relation 
alone does not fix the the entropy as a function of 
energy-momentum and particle number density. The 
reason is that the divergences are special derivatives 
and therefore the tensorial order of the coefficients a 
and Pa is lower than would be in case of a natural func- 
tional relation of the differentials. We will see that sev- 
eral different Gibbs relations are compatible with ((26)) . 
On the other hand, we may recognize, that (|^B)) is the 
balance of entropy constrained by the balances of the 
particle number density and the energy-momentum, 
where a and /3q are Lagrange multipliers^ 

The above formulas define the relation of kinetic 
theory and thermodynamics. Now we want to find the 
most general form of this correspondence. 



4 Thermodynamics and kinetic theory - 
natural frame 

According to the usual interpretation of the equilib- 
rium fields the absolute value of the vector (3a is the 
reciprocal temperature, its direction gives the veloc- 
ity field iia and a is the chemical potential over the 
temperature: 



T 



T 



(27) 



In this special case ()21l) is the classical equilibrium 
distribution of Jiittner [35] , the relativistic generaliza- 
tion of the non-relativistic Maxwell-Boltzmann distri- 
bution: 

h = e^^. (28) 

The momentum space integral defining the equilib- 
rium particle number density exist only if /3q is a time- 
like vector [33] . Therefore the particle number density 
four vector Nq in equilibrium is proportional to 13"^ 



(29) 



where A is a scalar field. 

This observation is the reason, why in local equi- 
librium one may expect zero diffusion current density 

= A^N'' = 0. In this case it is possible to define 
a velocity field by particle number density four vector 
assuming, that Nq = nu°' = Af]"", where n is the equi- 
librium local rest frame particle number density. In 
this case the vector is parallel to the (equilibrium) 
velocity field, and 



pa 



T 



1 



h:^Ay/^ (30) 



^ That later interpretation is important also beyond 
equilibrium, where H26p is an inequality, therefore a and 
Pa are Lagrange- Farkas multipliers [401141) . 
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are consistent definitions. We will call this choice of 
the velocity field as natural equilibrium frame. 

A similar calculation gives the equilibrium energy- 
momentum tensor. Together with the particle number 
balance they are written as 



(31) 



Tn"'' = / dt^pV/o = eu'^u^ - pA''\ (32) 

Here the diagonal form is the consequence of the 
special distribution function and we obtain the usual 
relations [22..39j : 



Hence the thermodynamic interpretation of the ki- 
netic quantities is complete and one may get the Gibbs- 
Duhem relation, by the total Legendre transformation 
of (123): 

dp = sdT + hdji. (41) 

Finally we should recognize, that according to ((36|) . 
the corresponding material is an ideal gas, providing 
the relation p = Tn. 

One can observe, that in the natural frame inter- 
pretation we have defined a convenient four velocity 
field and hence a rest frame of the material and all 
the thermodynamic quantities and relations are de- 
fined relative to that. 



h = A'Km^TK2 ( ^ I CT. 
p = fif, 



e = 3nT + mh 



K, (f) 



Substituting ^ and into (gS]) we get 



(33) 
(34) 

(35) 



(36) 



Therefore the direction of the entropy four vector 
is parallel to the natural velocity field. Substituting 
and dnil) into ^ results in 

(daS -f adafi - ^dae ] + 



daii'' ys + an- j^{e+p)j =Q. (37) 

We may require that the equality is valid for arbi- 
trary velocity fields and velocity gradients. Therefore 
the expression in parenthesis at the second line of the 
previous formula gives that the entropy density in the 
natural local rest frame is the function of the energy 
density and the particle number density s = s(e, h) 
and the partial derivatives are: 



ds 



ds 1 



dn t' de f 



(38) 



The consequent classical Gibbs relation for the differ- 
entials reads: 

de = Tds + (idfi. (39) 

The second line of ([57| gives the potential relation, fix- 
ing the Legendre transformation properties and defin- 
ing the pressure: 



p := Ts + fin — e. 



(40) 



5 Thermodynamics and kinetic theory - 
general frame 

We may assume, that the velocity field u° of an ideal 
fluid is not parallel to /3a- We will call this choice of the 
velocity field as general equilibrium frame. In this case 
/3a can be decomposed into parallel and perpendicular 
components: 

r ^ (42) 

where we have used the notation y := /3au"' and w"- := 

It is convenient to introduce also = -I- w". 
We can see that g°'Ua = 1 and w°'Ua = 0, moreover 
Waw"" < 1, because /3° and therefore g"" are timelike. 
Now the equilibrium distribution can be written as 



/o = e 



>J-P°("a+'JJa) 



(43) 



Then we can calculate the equilibrium particle num- 
ber four vector and also the energy momentum tensor, 
recognizing, that the relations of the different quanti- 
ties in natural and general frames are 



(44) 

where = —w°'Wa- 

Then the equilibrium particle number four vector 
A^Q and energy-momentum tensor Tq ^ can be written 
in the general frame as 

= nu" + nw" , (45) 

T^' = {e+p)g^g'-p5'^' 



q q 

e+p 



where 



Vl — w 
P = P, 

e + pw^ 



2 ' 



1 — ' 

g" = (e+p)u;^ 
E" = -I- (e-l-p)w". 



(46) 

(47) 
(48) 
(49) 

(50) 
(51) 
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The energy and particle number densities and the 
pressure can be calculated from the corresponding forms 
in the natural frame. Substituting ps)) and P5|) into 
(ESI) we get 

S- = ^^^—{nT-iin + E^g,), (52) 

where we have introduced := y/T^-w^ft. The en- 
tropy current density in the general frame is parallel 
to g"". The thermodynamic relations are calculated by 
substituting and (|55|) into initially assuming 
only that 5° = sg'': 

daSS + adaK - /JbS.T"" = 

— (r9as + tJ-daU - .gf,9a-B'' - pw''daUb) + 

^{Ts + tin^g,E''-p)=0. (53) 

The second line of the previous expression shows, 
that the entropy density in the natural local rest frame 
is the function of the particle number density, the 
energy-momentum density four vector and the veloc- 
ity s — s{E°' ^n,u°-'). The related partial derivatives 
are: 

ds _ _ /X ds _ ga_ _ds_ _ pwa_ , . 

The consequent Gibbs relation for the differentials is: 

g^dE"- = Tds + ^idn-pwadu". (55) 

According to this relations, the change of the velocity 
directly contributes to the thermodynamic quantities. 

The expression in the parenthesis of the third line 
of ([5^ is the potential relation, fixing the Legendre 
transformation properties as: 

p^Ts + ^in- gaE"-. (56) 

Then the Gibbs-Duhem relation follows: 

dp + E'^dga = sdT + ndp, - pUadw" . (57) 

Finally ((52|) shows, that the corresponding mate- 
rial is an ideal gas, providing the relation p = Tn. 
This equation of state does not change, is the same 
both in the natural and general frames. 

One can check the above expressions directly, rec- 
ognizing, that 

s{E'',n)^-s{ng,E-ga), (58) 
9 

p{g-,T,f,)^p(-,l^]. (59) 

Here g = \J\ — v? = \fg^ ■ However, these relations 
alone, without the constrained entropy balance (j26p . 
are not enough to fix a unique form of the differentials. 



One may see more the role of the pressure and the 
velocity change - the acceleration - related term in the 
Gibbs relation (|55p , if the energy-momentum density 
vector is written in the form i?" ~ eu'^ + q°'. In this 
case we get 

Tds + fidn — gadE°' + pwadu'^ = 

{ua + Wa)d{eu°- + q"-) + pwadu°- = 

de + Wadq" + ((e -|- p)wa - qa)du'' . (60) 

Therefore taking into account (jSOp we can free our- 
selves from the direct pressure and velocity depen- 
dence and get 

de -t- Wadq"' = Tds + p,dn. (61) 

In this relation the velocity field dependence is not 
apparent. 



6 Summary and conclusions 

It is well known, that thermodynamic relations are 
fixed to the material and therefore they are best ex- 
pressed by rest frame quantities. This is a hidden as- 
pect in nonrelativistic hydrodynamics |44| , but an ap- 
parent property of relativistic theories. In thermody- 
namic equilibrium there is a single, distinguished ve- 
locity field, and a corresponding rest frame, where all 
motion related physical quantities are simple. 

In this paper we have investigated whether a novel 
concept of relativistic thermodynamics based on a gen- 
eralized Gibbs relation is compatible with kinetic the- 
ory. Kinetic theory distinguishes a natural frame, a ve- 
locity field in equilibrium, fixed by the direction of 13a- 
The entropy and particle number density four vectors 
are parallel to this velocity and the energy-momentum 
density generated accordingly. Therefore this natural 
frame corresponds both to Eckart and Landau-Lifshitz 
frames of dissipative fluids ^SIHS] . In a natural frame 
the thermodynamic relations are very similar to the 
nonrelativistic ones. 

However, one may introduce a general frame, a ve- 
locity field independently of /3a, and investigate the 
consequent thermodynamic relations. The calculations 
are transformations of the equations calculated in the 
natural frame, but the resulted thermodynamic rela- 
tions contain additional terms. In particular the energy- 
momentum vector appears as a natural variable - sim- 
ilarly to (0) - and the velocity field plays an explicit 
role in the relations. 

Therefore we have shown that kinetic theory is 
compatible with the generalized Gibbs relation sug- 
gested to resolve stability problems of hydrodynamics 
and paradoxes related to the temperature of moving 
bodies. The direct velocity dependence of the Gibbs 
relation (j55p in the general frame generalizes our pre- 
vious results, where we either assumed that the en- 
tropy production is independent of acceleration |24j 
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or considered a thermodynamic body moving with a 
uniform velocity T. 

If ([55)1 (and ([55]) ) are introduced in a hydrody- 
namic theory, the concept of local equilibrium is mod- 
ified. Because the constitutive relations of dissipative 
fluids (viscosity, heat conduction, etc..) express the 
tendency toward local equilibrium, the entropy pro- 
duction and thermodynamic fluxes and forces can be 
modified, too. 
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